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The effect of pressure on the superconducting transition temperature (Tc) of the Ca-intercalated 
graphite compound CaCe has been investigated up to ~ 16 kbar. Tc is found to increase under 
pressure with a large relative ratio /STc/Tc of ~ +0.4%/kbar. Using first-principles calculations, we 
show that the positive effect of pressure on Tc can be explained within the scope of electron-phonon 
theory due to the presence of a soft phonon branch associated to in-plane vibrations of the Ca atoms. 
Implications of the present findings on the current debate about the superconducting mechanism in 
graphite intercalation compounds are discussed. 



PACS numbers: 74.70.Ad, 74.62.Fj, 74.25. Kc, 74.62.-c 



I. INTRODUCTION 

The discovery of superconductivity in Yb- and Ca- 
intercalated graphitesii^, with significantly higher T^s 
than found previously for alkali-metal graphite interca- 
lation compounds (GICs) re-initiated the debate on the 
long-standing puzzle of the origin of superconductivity in 
GICs'^. Possible pairing mechanisms under consideration 
range from unconventional exciton- or plasmon-mediated 
pairing"* to the more conventional electron-phonon {e-ph) 
coupling mechanism^i^. 

The core of the debate on the origin of supercon- 
ductivity in the GICs are the so-called interlayer elec- 
tronic bands which cross the Fermi level (Ep) besides 
the graphite n* bands. Since the coupling of the tt* band 
electrons to the graphite phonons is not sufhcient to yield 
the observed high Tc's^'^, the role of the interlayer bands 
and their degree of filling is considered to be essential 
for superconductivity with increased Tc's^. Csanyi et al. 
proposed that the electrons in the interlayer bands can 
be considered as nearly-free 2-dimensional (2D) electrons 
propagating essentially between the graphene layers^, 
and pairing can be mediated by acoustic plasmons in 
these 2D metallic slabs^. However, there is a growing 
body of evidence that the interlayer bands have suffi- 
ciently strong coupling with both in-plane intercalant 
and out-of-plane graphite phonon modes to allow for the 
relatively high Tc's within the scope of a standard e-ph 
coupling mechanism^. This scenario gained some exper- 
imental support from recent penetration depth'^ , specific 
heat^", tunneling*^ and Ca isotope effeclji^ experiments 
in CaCg. 

Interestingly, the latter experiments showed a Ca iso- 
tope effect with aca = 0.5, suggesting a dominant role of 
Ca-related modes in the superconducting pairing^^. This 
finding, however, is at odds with the results of ab-initio 
calculations^ in which isotope exponents aca ~ etc ~ 
0.25 have been concluded. The fact that the out-of-plane 
phonon modes of the host layers cannot be ignored gains 
substantial support from recent investigations on other 



layered compounds, e.g. CaAlSi (Tc ^ 8 K), where Ca 
atoms are located in between Al-Si honeycomb sheets 
(AIB2 structure type)^'^. Hence, a deeper understanding 
of the superconducting pairing, and of the nature and 
role of the electronic and vibrational states involved, will 
provide useful insight into superconductivity in the GICs 
as well as other recently discovered " high- Tc" supercon- 
ductors. 



The investigation of the pressure dependence of Tc is a 
key experiment to test and compare these theoretical pre- 
dictions. The degree of filling of the interlayer bands, cru- 
cial for the acoustic plasmon pairing mechanism, depends 
not only on the charge transfer from the intercalant, but 
also on the separation of the graphene sheets^. By apply- 
ing hydrostatic pressure the graphite layer spacing, and 
hence the energy of the interlayer bands, can be continu- 
ously tuned without affecting the chemical composition. 
On the other hand, a comparison of the pressure depen- 
dence of Tc with the results of first-principles calculations 
of the electronic and vibrational properties of the GICs 
could support or rule out the hypothesis of an e-ph me- 
diated coupling mechanisroi^. 



In this paper, we report the effect of pressure (P < 16 
kbar) on Tc of CaCe and compare our results with ab- 
initio calculations of the electronic and vibrational prop- 
erties carried out as a function of pressure. Based on our 
calculations, we argue that the increase of Tc with pres- 
sure can be understood within the scope of e-ph coupling 
as a consequence of a softening of phonon modes involv- 
ing Ca in-plane vibrations, while it is at odds with an 
acoustic plasmon mechanism. We also show that there is 
some discrepancy between the Tc measured from experi- 
ment and that predicted by isotropic, harmonic Migdal- 
Eliashberg theory, and discuss possible ways to improve 
the agreement between experiments and theory. 
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II. EXPERIMENT 

CaCg samples were prepared as described in detail 
elsewhere^iiiS. The Tc's of three samples were determined 
from the temperature dependence of the dc (samples SI 
and S2) and of the ac magnetic susceptibility (sample 
S3). The variation of Tc between our samples is less than 
~ 0.1 K. All samples show a sharp superconducting tran- 
sition with a width ATc - 0.15 K between 10% and 90 % 
of the diamagnetic signal, indicating good sample qual- 
ity. Cu-Be piston-anvil-type pressure cells were used to 
apply quasi-hydrostatic pressures with silicon oil as the 
pressure transmitting medium. To monitor the pressure, 
Tc of Pb (99.9999%) was determined simultaneously. 

At ambient pressure, the onset of the superconducting 
transition of our samples is ^ 11.4 K, consistent with pre- 
vious reports^.. The superconducting transition is clearly 
shifted to higher temperatures as the pressure is raised 
(see Fig.[lja)). Up to 16 kbar there is no indication of an 
abrupt change of Tc or a narrowing or broadening of the 
transition, as frequently observed in other superconduct- 
ing GICs, such as e.g. KHgC4, KCs, RbCs, and ascribed 
to a pressure-induced change of the staging or of the in- 
tercalant sublattice^^. No hysteresis in Tc as well as in 
the shape of x{T) was observed between the data taken 
with increasing and decreasing pressure. Thus the possi- 
ble shear stresses due to solidification of pressure medium 
do not affect significantly the superconducting properties 
of the samples. 

The pressure dependence of Tc for three samples is 
shown up to 16 kbar in Fig. [ijb). For all samples, Tc 
increases under pressure almost linearly with a slope of 
0.042 - 0.048 K/kbar. Depending on the sample, Tc and 
the ratio dTc/dP vary slightly, but the relative change 
{l/Tc)dTc/dP is the same within error bars for all sam- 
ples and amounts to « -1-0.4 %/kbar^^. Such a pres- 
sure dependence of Tc is adverse to the behavior of many 
sp superconductors in which often a negative pressure 
dependence of Tc is found. Also it is in contrast to 
other intercalated graphites, such as KHgCs, RbHgCs 
and KTI1.5C4 whose Tc decreases under pressure with a 
rate of -2.2 - -3.8 %/kbari^. 



III. THEORY 

To gain more insight into the possible coupling mech- 
anism we performed a detailed Density Functional Per- 
turbation Theory (DFPT)''^^ study of the electronic and 
vibrational properties of CaCg and explored the pres- 
sure range up to 500 kbar to investigate whether lattice 
instabilities occur. In our calculations we employed ul- 
trasoft pseudopotentials^ with a generalized Gradient 
Approximationi^ for the exchange-correlation functional. 
The eigenf unctions were expanded on a plane-wave basis 
set3 with a cut-off energy of 30 Ryd and 300 Ryd for 
the wavefunctions and the charge density, respectively. 
We used a (8)'^ Monkhorst-Pack grid and a 0.06 Ryd 
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FIG. 1: (Color online) (a) Temperature dependence of the 
susceptibility for CaCs (SI) at different pressures. The num- 
bers next to the data and in the bracket corresponds to the 
applied pressure (kbar) and the sequential order of the mea- 
surement runs. Tc is determined as the temperature where 
the extrapolation of the steepest slope of x{T) intersects the 
extrapolation of the normal state x{T) to lower temperatures, 
(b) Pressure dependence of Tc for three CaCe samples. The 
filled and open symbols are data taken at increasing and de- 
creasing pressure, respectively. The dashed lines are guide-to- 
eyes. The inset shows the relative change of Tc with pressure. 



Methfessel-Paxton smearing, which led to convergence 
of better than 0.1 mRyd for the total and 5 cm~^ for 
the F-point frequencies. Phonon dispersion curves were 
obtained by a Fourier interpolation of the dynamical ma- 
trices calculated on a 4^^ grid in q-space for P=0, 50 and 
100 kbar and a 2^ grid elsewhere. 

The electronic structure and phonon dispersion at am- 
bient pressure are in excellent agreement with previous 
results^. The Equation of State calculations yielded equi- 
librium lattice constants {arh = 5.16 a.u., rhombohedral 
angle = 49.90°) in very good agreement with the ex- 
perimental data {arh — 5.17 a.u., 6 — 49.55°), and a 
bulk modulus {Bq « 103 GPa) larger than that of pure 
graphite, but close to that of other GICs^^. The com- 
pressibility remains anisotropic {kc/ka ~ 9) but is signif- 
icantly reduced as compared to non-intercalated graphite 
[see Fig. 2(b)]. 

First of all, we notice that, in general, the coupling 
strength between electrons and any bosonic modes can be 
approximated by A cx N{Ef)(I^) / {to'^) , where N{Ef) is 
the electronic density of states at Ep, (J^) is the electron- 
boson matrix element, and (w^) is the square averaged 
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FIG. 2; (Color online) Ab-initio pressure dependence of se- 
lected structural, electronic and vibrational properties of 
CaCe- (a) Frequencies of the F- (open symbols) and X- (solid 
symbols) point phonons. Imaginary frequencies are shown as 
negative. The labels indicate the atom giving the dominant 
contribution to the phonon eigenvector, (b) Variation of the 
in-and out-of-plane lattice constants and of the relative ratio, 
(c) N{Ef) and position of the bottom of the interlayer bands 
with respect to Ef- 



frequency for the relevant bosonic modes. Our calcula- 
tions show that N{Ep) decreases with pressure because 
and the interlayer band is emptied and the tt* bandwidth 
increases, (cfr. Fig. [2^c)). In the following, we therefore 
investigate the alternative option to enhance A (and Tc) 
by a reduction of the frequency of the relevant bosonic 
modes. 

First, we consider electron-phonon coupling. Although 
unusual, the softening of a particular phonon mode, usu- 
ally an indication of an incipient lattice instability, can in- 
duce a pressure increase in T^. The calculated phonon 
modes for CaCe behave differently as a function of pres- 
sure depending on their eigenvector as seen, 6.17. in the 
case of the T-point phonons (see Fig. O^a)): whereas the 
C in-plane modes (w > 100 meV) and the Ca out-of- 
plane mode (w ~ 30 meV) harden with pressure, the C 
out-of-plane mode (50 < uj < 100 meV) are almost unaf- 
fected. On the other hand, the lowest-lying optical mode 
at r, and an acoustical mode at X (Fig. 2]), both mainly 
involving in-plane Ca vibrations, are considerably soft- 
ened with pressure. At higher pressures {P > 120 kbar), 
these modes drive the system unstable. Frozen-phonon 
calculations reveal that this mode couples to both the in- 
terlayer and the w* bands, and that the mechanism which 
causes its instability is similar to that giving rise to the 
anharmonicity of the mode in MgBz^l. 

To estimate the pressure dependence of the partial Aca 
associated with the Ca in-plane vibrations, we employ a 
simple approximation based on the Hopfield formula: Xi 
= N{EF){lf)/M,[u;f), where {If) is the e-ph matrix ele- 
ment. Mi is the phonon mass, and (uf) is square averaged 
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FIG. 3: (Color online) Cross section of the Fermi surface of 
CaCe at three different pressures (P = 0,100,250 kbar), in 
a plane orthogonal (a) and parallel (b) to the c-axis. The 
spherical, 3D Fermi surface centered around F is due to the 
interlayer band, while the 2D tt* bands form distorted cylinder 
along the c axis. The position of the special k points in this 
plot are relative to the P = lattice constants. 



phonon frequency. Assuming that (1^^) constant and 



can be approximated by the square of the lowest- 
lying optical phonon frequency at F, we find a significant 
increase of Aca by 20% at P = 50 kbar, and 60% at P 
— 100 kbar with respect to its zero-pressure value. Since 
the low energy phonon modes give dominant contribu- 
tion to the coupling^, the increase of Aca can be strong 
enough to overcome the reduction of N{Ep) and as a 
result, increase Tc. 

In contrast, it appears that the other proposed pairing 
mechanism, namely acoustic plasmon-mediated pairing 
in the interlayer bands^i?-^ can hardly be reconciled with 
the observed pressure dependence of Tc- The basic as- 
sumption of this model is the existence of 2D metallic 
layers sandwiched by dielectric layers. If this is the case, 
due to the significant anisotropy of the Fermi velocity 
and incomplete screening, the plasmon dispersion devel- 
ops a low-energy acoustic branch, Qpi{c{) ^ q, which can 
provide additional pairing routes for superconductivity^. 

However, the electronic band structure (cfr. Ref. (|) and 
the corresponding Fermi surface (FS) plot (see Fig. [3]) 
reveal that the interlayer bands are 3D in character. Ad- 
ditionally, there is a substantial charge transfer to the 
graphite tt* states. Both findings question the basic as- 
sumptions of this model already at ambient pressure. 
The anisotropy of the Fermi velocity emerges almost ex- 
clusively from the graphite n* bands, which display a 
warped cylindrical FS parallel to the c axis. With pres- 
sure, the dispersion of the tt* bands along the c axis in- 
creases and the 3D character of the electronic structure 
grows. When the system becomes more 3D, a gap in the 
acoustic plasmon dispersion is developed^, and the con- 
ventional plasmon dispersion, npi{q)=flpi-\-0{q^) {^Ipi 3> 
ksTc) is fully recovered. Thus, even if any coupling due 
to acoustic plasmons is assumed, its strength will de- 
crease with pressure, adverse to our experimental obser- 
vation. 

Finally, we check whether the increase of A due to 
the phonon softening is compatible with the observed in- 
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where <uji„> is the logarithmic averaged phonon fre- 
quency and /i* is the Coulomb pseudopotential. Insert- 
ing the calculated <cji„>, (305 K at P = 0, 300 K at 
P = 50 kbar, and 230 K at P = 100 kbar), and setting 
fi* — 0.145, which reproduces the experimental Tc at zero 
pressure, we estimate that Tc should increase from 11.4 K 
to 13.5 (12) K from P = to P = 100(50) kbar. There- 
fore the effect of the phonon softening for the in-plane 
Ca vibration is strong enough to overcome those of the 
decrease of DOS and <w;„>, and as a result, increases 
with pressure. We also predict that the phonon soften- 
ing will drive the system unstable for P > 120 kbar, thus 
imposing a theoretical limit to the maximum attainable 
Tc. 



IV. DISCUSSIONS 



FIG. 4: (Color online) (a) Phonon frequencies and (b) density 
of states of CaCe for P = (black), P=50 (green), P = 100 
kbar (orange), along selected directions in the rhombohedral 
unit cell; the line F — X is contained in the graphene planes, 
while L — F is orthogonal to it. (c) Eliashberg function a^F{Q) 
and frequency-dependent electron-phonon coupling A(a;) — 
2j^a'^F{Q,)/ndQ. 



crease of Tc under pressure by calculating the e-ph spec- 
tral function a^F^uj) and the frequency-dependent e-ph 
coupling parameter X{uj) = 2 /p a^F{U)il,~^dn. Instead 
of the maximum experimental pressure (P w 16 kbar), 
we choose to run these calculations at two higher pres- 
sures, P = 50 and 100 kbar. We notice, in fact, that the 
experimental variation of Tc up to P w 16 kbar is too 
small to be resolved by the accuracy of DFT calculations. 
Therefore, we choose two pressures where the calculated 
frequencies are all real {i.e. the system is structurally 
stable), and where our approximated Hopfield equation 
predicts a sizable increase in A; the second point (P — 
100 kbar) is chosen closed to the predicted structural 
transition, in order to get an estimate of the maximum 
attainable Tc in this system. 

The e-ph spectral function a^F(u;) displays three 
peaks well separated in energy, corresponding to differ- 
ent phonon eigenvectors. The largest contribution to the 
total e-ph coupling comes from the in-plane Ca phonons, 
which increases Aca = 0.42 to Aca = 0.62(0.44) at P = 
100 (50) kbar. The second largest contribution from the 
C out-of-plane vibrations remains essentially unchanged, 
since the modes at ~ 70 meV are insensitive to the pres- 
sure. The total A increases from 0.84 at P = to 0.86 
at P = 50 kbar and 1.03 at P = 100 kbar, which results 
entirely from the soft Ca in-plane modes. 

We estimate the effect of pressure on Tc using the 
McMillan formulaSi: 



1.2 



■ exp 



-1.04(H- A) 



After finishing the present work, an experimental study 
up to P w 160 kbar was posted on the cond-mat 
archive^^. Here a transition of Tc at ~ 80 kbar was ob- 
served, accompanied by lattice softening, thus confirming 
our theoretical prediction of an upper bound on Tc due 
to lattice instability. In Fig. [5] we plot the experimen- 
tally measured Tc, from our work and Ref. [2^, together 
with the calculated Tc as a function of pressure. We no- 
tice that, even though the positive pressure dependence 
is well reproduced by our ab-initio calculations, there are 
still some discrepancies between experiment and theory, 
which we shall now discuss in detail. 

First of all, while experimentally the (structural) tran- 
sition is found at 80 kbar, we predicted it to take place 
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FIG. 5; (Color online) The calculated Tc (open squares) as 
a function of pressure together with the measured Tc in this 
work (same as Fig. [2]) as well as Ref. 25 (green triangles) . The 
dashed line indicates the point at which ab-initio calculations 
predict a dynamical instability of the lowest X-point phonon. 
The inset shows a pressure dependence of the calculated e- 
ph coupling constant, A and the logarithmic averaged phonon 
frequency, <tJ!„>. Both quantities have a strong non-linear 
behavior. 
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at 120 kbar. It is well known that the structural tran- 
sition pressure for the GICs is very sensitive to several 
experimental details. For example, the previous pressure 
experiments on Ist-stage KCg showed large hysteresis of 
the transition pressure as well as significant time depen- 
dence of the order of days or weeks^^. In particular, we 
note that in the GICs Tc can show a strong anomaly even 
before a structural transition takes place. For example, 
in KCg, an intercalant structural transition occurs at ~ 
14 kbar—, but a sudden increase of is observed already 
at ~ 4 kbar—. On the other hand, this kind of error is 
also within the accuracy of DFT. While it could possi- 
bly be reduced by a more careful convergence of the X 
point frequency with respect to k points sampling, there 
would always be an uncertainty of 10 - 20 kbar related 
to the choice of the exchange and correlation functional 
or of the basis selj^. However, even if DFT can fail to 
predict the exact transition pressure, it is usually much 
more accurate in describing structural trends as a func- 
tion of pressure. Therefore, it is quite surprising to notice 
that the behavior of Tc with pressure predicted by theory 
deviates from that measured experimentally: in particu- 
lar, theory predicts a much slower increase in Tc (~ 0.02 
K/kbar) than experiment, and a slightly non-linear be- 
havior, with a stronger increase in Tc at higher pressures. 

Such discrepancies indicate that the harmonic, 
isotropic e-ph coupling theory may not be sufhcient to 
describe quantitatively the behavior of CaCg and other 
GICs. This confirms what was shown by a very recent, 
puzzling experiments on the isotope effectlS,, which re- 
ported surprisingly high isotope exponent for Ca, a(Ca) 
« 0.5, much higher than the value a(Ca) = 0.24 predicted 
by theory^. Assuming a nontrivial contribution from the 
C phonon modes, the total isotope exponent will exceed 
the BCS value of 0.5. Also, discrepancies between the 
isotropic, harmonic Migdal-Eliashberg theory and exper- 
iment are found in the specific heat and the upper critical 
fields'^'^. Further theoretical and experimental works are 
highly desirable along this direction. 

The first possibility to improve the agreement between 
theory and experiment would be taking into account an- 
harmonic effects or nonlinear coupling of the in-plane Ca 
phonon modes32-. In graphite intercalation compounds, 
the vibrational modes associated to the intercalant are 
extremely soft, and transitions to different in-plane sub- 
lattice as well as staging are frequently observe d 
which indicates that these modes are strongly anhar- 
monic. 

The second approximation which may not be justified, 
in the present case, is that of an isotropic electron-phonon 
coupling. As we have previously discussed, the Fermi 



surface of CaCg is highly anisotropic and is formed by 
two bands, which have very different electronic origin 
(tt* and interlayer), and are hence coupled to different 
phonon modes. Therefore, different sheets of the Fermi 
surface would give very different contributions to the to- 
tal electron-phonon coupling, possibly giving rise to a 
"smeared" multi-gap superconductivity^. 

Finally, we would like to compare our results with a 
very recent ab-initio calculations-^. Here the authors 
study, using the same codo22., but different pseudopo- 
tentials, the behavior of CaCe as a function of pressure. 
Similarly to us, they find a structural transition driven 
from the softening of the X-point phonon, at some pres- 
sure between 70 and 100 kbar, thus lower than ours and 
close to the experimental value. Between and 50 kbar, 
they find an increase in A of only 0.015, and estimate 
an increase in Tc of ^ 0.3 K; they do not calculate the 
Tc at any other pressure besides and 50 kbar, so it is 
not possible to compare directly the pressure behavior of 
the two sets of calculations, but the overall qualitative 
picture seems to be the same; the numerical differences 
in the transition pressure and A are probably due to the 
use of different pseudopotentials and k-point sampling. 
To explain the difference between theory and experiment, 
the authors of this paper propose still another possibil- 
ity, that a continuous staging transition takes place under 
pressure. This is of course a possibility which should be 
taken into account, but a definitive answer in this sense 
could come only from X-ray or neutron diffraction mea- 
surements under pressure. 

In conclusion, we presented measurements and ab- 
initio calculations of the pressure dependence of Tc in 
CaCg. We demonstrated that the positive pressure de- 
pendence of Tc can be understood within an e-ph sce- 
nario due to softening of the in-plane Ca phonon modes, 
while it appears not to reconcile with the acoustic plas- 
mon mechanism. In view of new experimental findings, 
we also discuss which effects beyond the isotropic, har- 
monic Migdal-Eliashberg theory would have to be taken 
into account to obtain a full quantitative and qualitative 
agreement between theory and experiments in CaCe, and 
possibly other graphite intercalation compounds. 
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